Current constraints on the cosmic growth history 
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We present constraints on the cosmic growth history with recent cosmological data, allowing for 
deviations from ACDM as might arise if cosmic acceleration is due to modifications to GR or in- 
homogeneous dark energy. We combine measures of the cosmic expansion history, from Type la 
supernovae, baryon acoustic oscillations and the CMB, with constraints on the growth of structure 
from recent galaxy, CMB and weak lensing surveys along with ISW-galaxy cross-correlations. De- 
viations from ACDM are parameterized by phenomenological modifications to the Poisson equation 
and the relationship between the two Newtonian potentials. We find modifications that are present 
at the time the CMB is formed are tightly constrained through their impact on the well-measured 
CMB acoustic peaks. By contrast, constraints on late-time modifications to the growth history, 
as might arise if modifications are related to the onset of cosmic acceleration, are far weaker, but 
remain consistent with ACDM at the 95% confidence level. For these late-time modifications we 
find that difi'erences in the evolution on large and small scales could provide an interesting signa- 
ture by which to search for modified growth histories with future wide angular coverage, large scale 
structure surveys. 



I. INTRODUCTION 

A central challenge facing cosmology today is to ascer- 
tain the origin of current accelerated cosmic expansion 
(see for example [IHl])- Precision measurements of geo- 
metric distances to astrophysical objects - the luminosity 
distance to Type la supernovae (SNla), the ratio of an- 
gular diameter and radial distances from baryon acous- 
tic oscillations in the 3D galaxy distribution (BAO) and 
the angular diameter distance to last scattering of the 
Cosmic Microwave Background (CMB) - have provided 
consistent evidence that the universe underwent a tran- 
sition to acceleration rather recently, when the universe 
was roughly half the size it is today. 

Distance measurements, however, purely probe the ho- 
mogeneous expansion history. This is only part of the 
story if acceleration derives from a modification to Gen- 
eral Relativity (GR) on large scales, or a new type of 
matter, rather than a cosmological constant (A). If this 
diverse range of possible causes for acceleration are con- 
sidered, a clearer discrimination between them is gained 
from also considering the growth of density fluctuations, 
and the large scale structures they seed. 

Recently there has been significant interest in the po- 
tential of upcoming surveys to constrain the growth of 
structure including multi-frequency imaging, galaxy and 
weak lensing measurements |4|-|l8{ , spectroscopic surveys 
to measure the peculiar velocity distribution jH, [13, lisl - 
[2^ . and 21cm intensity surveys H^. A principle focus 
is the application auto- and cross-correlations of galaxy, 
weak lensing and velocity fields to ascertain if signatures 
of deviations from a standard scenario, of GR and a ho- 
mogeneous dark energy component, can be detected in 
the growth history. 

Existing measurements of the growth of structure can 
already provide constraints on deviations from this stan- 
dard scenario, however. Studies have included consider- 
ing phenomenological models of modified gravity l24-|32 
the presence of dark energy inhomogeneities |33| - |40| | and 



specific modified gravity theories such as the DGP |4l| 
scenario and /(i?) gravity [H-il. 

In this paper we provide a current reference point, of 
how well existing cosmological data can constrain the 
growth of structure, against which the future survey tar- 
gets can be framed. We use a combination of both the 
current cosmological distance measures and auto- and 
cross- correlations of large scale structure measurements 
to place constraints on the growth history when it is not 
required to be dictated by GR. Type la supernovae from 
the Supernovae Cosmology Project [H^l and the CfA3 
sample [55|, the joint analysis of BAO acoustic signa- 
tures in the Two-Degree Field (2dF) and DR7 Sloan 
Digital Sky (SDSS) surveys [s^, and the WMAP 7-year 
(WMAP7) CMB temperature and polarization spectra 
[57| provide constraints on the cosmic expansion history. 
The WMAP 7- year CMB data, in combination with the 
matter power spectrum from the SDSS DR7 release [sl] , 
Integated Sachs Wolfe (ISW)-galaxy cross-correlations 
with 2MASS and SDSS Luminous Red Galaxy (LRG) 
galaxy surveys f59| and the COSMOS weak lensing sur- 
vey |60[ constrain the growth history. 

We employ a phenomenological framework to describe 
how the growth history might be modified, rather than 
choosing to investigate constraints on specific modified 
gravity models currently proposed in the literature. A 
constant or simply evolving equation of state is used to 
investigate observational signatures of the origin of cos- 
mic acceleration in the background expansion history. 
Here we use simply time- and scale-varying models to 
bridge the link between observation and theory for the 
growth history. 

In section |lT] we discuss phenomenological changes to 
the standard growth history, modifying the perturbed 
Einstein equations through the Poisson's equation and 
the relationship between the two Newtonian potentials. 
We also consider four model-independent parameters 
commonly used to describe cosmic growth history. In 
section Hill we discuss the data sets used in the analysis 
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and in section IIVI summarize the effects that the modi- 
fied growth history can have on the observed large scale 
structure correlations. In section |Vl we summarize the 
results of Monte Carlo Markov Chain (MCMC) analyses 
to place constraints on the cosmic growth history, and 
draw together our findings and discuss the implications 
for future work in section IVTl 



where = Si + 3H{l + w)9i/k'^ is the rest-frame density 
perturbation of matter species i,'H{a) = a/ais the Hub- 
ble expansion factor, where dots represent derivatives 
with respect to conformal time and Ui the anisotropic 
shear stress. Equation [5] shows that the two potentials 
become effectively equal in the matter and dark energy 
dominated eras, when there are negligible anisotropic 
shear stresses. 



II. MODELING THE GROWTH OF 
STRUCTURE 



B. Deviations from the standard growth scenario 



We consider the conformal Newtonian gauge to de- 
scribe the metric 



-a(r)2 [1 + 2iA(x, t)]dT^ + a(r)2 [1 - 20(x, t)]dyi^ (1) 



We parameterize deviations to the growth history by 
modifying equations ([4]) and (O with two scale and time- 
dependent functions, Q{k,a) and R{k,a), 



where a is the scale factor (a = 1 today), r is confor- 
mal time, X are comoving coordinates, and tp and (j) are 
the two Newtonian potentials describing temporal and 
spatial perturbations to the metric respectively. We as- 
sume units in which the speed of light c = 1. The 
spatial element can be broken down into a radial, line 
of sight, component, written in terms of the comoving 
distance, x = fa '^'^i ^^^"^ ^ 2D angular element dfi^, 
dx^ = (dx^ + r-(x)^df^^) , where r(x) is the comoving 
angular diameter distance. 



A. Growth in General Relativity 

The growth of inhomogeneities is characterized by 
the evolution of the fractional over-density distribution, 
S{k,a) = p{k,a)/p{a)~l, and the divergence of the pecu- 
liar velocity, 6{k,a), where k is the comoving wavenum- 
ber. 

We employ energy-momentum conservation in the per- 
turbed fluid to obtain the fluid equations for each type of 
matter, with the density and peculiar velocities evolving 
according to [Slj . 

5 = -{l+w){e-i^)-m{cl-w)5 (2) 
6 , , d w 9 c2 ^ 

B ^ ' k^ l + wB {l + w) ^ ' 

Here w = P/p is the matter's equation of state and = 
dP/dp the sound speed. 

To wholly describe the growth history of density and 
velocity perturbations (up to initial conditions) two fur- 
ther equations arc required. We use the Poisson equation, 
combining the time-time and time-space components of 
the perturbed Einstein equations, and the anisotropic 
space-space component, 

k^(f> = -AirGa^^piA, (4) 
= -127rGa2^p,(l + u;,)^. (5) 



i 

^-R(j) = -l2nGQa^Y.p,{l + w)j. 



(6) 



(7) 



We assume throughout that matter remains minimally 
coupled to gravity so that the fluid equations ^ and ([3]) 
remain unchanged. 

In the context of modified gravity scenarios, equation 
dS]) describes a modified Poisson equation in which the 
gravitational potential responds differently to the pres- 
ence of matter, while ([7]) allows an inequality between 
the two gravitational potentials, even at late times when 
anisotropic shear stresses are negligible, with tp « Rcj). 
Note that the time-evolving modifications (|6]) and ([7]) 
are consistent with the predicted evolution of the super- 
horizon curvature fluctuations, as discussed for example 
in [13, © simply provides the constraint equation by 
which both (f> and ip superhorizon solutions can be found. 

Appendix lAl summarizes how the modified growth his- 
tory is incorporated into the publicly available CAMB 
code [13 , used to predict the effect of the model on CMB 
and galaxy statistics. 

The naming conventions for these phenomenological 
modifications have been somewhat varied in the litera- 
ture: Q has been defined in other works as Q [gl], / [sst . 
Geff/G [13, Hi, 9 H and /xJllLP, [HjWhile i? is 
defined as IMI -h r?) [6| , I/t? [1111371/7 [ill, 0, and 
1 + 137, [23, do, [S^- We choose Q and R here in order 
to avoid confusion with the synchronous gauge metric 
variable, 77, discussed in appendix [A1 and growth history 
variables, go, /, and 7, introduced in section Hi CI 

Theoretical models of modified gravity can be de- 
scribed using ^ and ([7]), for example exhibiting time- 
dependent, but largely scale-independent variations, as 
with DGP on sub-horizon scales, in the quasi-static 
regime, [68l-l70|. or with both time and scale dependence, 
for example in f{R) theories [zO, [zH- We consider two 
phenomenological models in our analysis that allow for 
both time and scale dependence in the modifications. 

The first (Model I) parameterizes Q and R with mono- 
tonically evolving, time- and space-dependent functions 
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of a and k, with 
Q{k,a)-l = 



he 



i?(fc,a)-l = [i?oe"''/''<^ +i?oo(l - e^'^/'^'^) - ij a'*(8) 

where Qq,Ro and Qoc,Roo are, respectively, the asymp- 
totic values of the modification on superhorizon and sub- 
horizon scales today. In specific theories, Q and R arc 
typically a monotonically varying function of time and 
scale. Here, for simplicity, we assume a constant power 
law index, s, to parameterize the time-variation of the 
modification, kc is a comoving transition scale between 
the large and small scale behavior that allows us to de- 
couple observations on large and small scales to assess 
which are providing the principal constraints. 

We impose two theoretical priors on the model, one 
that Q{k,a) > 0, so that matter is attracted into over- 
densities, and R{k,a) > — 1 to ensure that null geodesic 
paths are bent inwards by gravitational potentials. 

Modifications to gravity could be mimicked by modi- 
fications to the growth of inhomogeneities in the 'dark 
sector' [tI, [t^, since the effect of any additional in- 
homogeneities is only inferred indirectly through their 
gravitational effect. Such modifications could include 
the presence of dark energy density, peculiar velocity 
and anisotropic shear stress fluctuations or non-minimal 
matter-dark energy interactions that contribute to the 
right hand side of Einstein's equations. In this context, 
the inhomogeneities would give rise to effective Q and R 
parameters. 



Q 



R. 



■eff 



= 1 



= 1 



de 



Spde 



(9) 



(1 -I- Wde)crde + (-S^) (1 + Wm)(Jr, 



Pm^m + Pde^d 



(10) 



where pm^m is the net inhomogeneity from Standard 
Model and cold dark matter, pde^de describes the ad- 
ditional dark energy and coupled dark energy-dark mat- 
ter inhomogeneities, and ade is a dark energy anisotropic 
shear stress. 

If the dark component producing the inhomogeneities 
is non-relativitic, then it is extremely difficult to sustain a 
substantial anisotropic shear stress (see for example (t^). 
If, in addition, at early times, the dark energy density in- 
homogeneities are small in comparison to those of normal 
matter, then one may reasonably assume Reff = 1- 

The second model (Model II) fixes R = 1 and modifies 
the growth history purely through including an extra, 
additive, time- and space-dependent contribution to the 
Poisson equation 



Pn 



-fc/fc. 



where A^h parameterizes the amplitude of the modifi- 
cation on horizon scales (fc = H^^) today, and rix 



describe the temporal and large scale spectral evolution, 
and on scales k > k^ the modification is suppressed (fc^ 
would relate to an effective comoving sound horizon, if 
the modification was due to an additional component). 
As an illustration, a modification with = ns/2 {hs 
is the inflationary spectral index), Sx = and kx = oo 
would mimic large scale CDM density perturbations in 
the CDM dominated era. 

In the presence of pressureless matter this modification 
is equivalent to 



= 1 



-k/k^ 



ilmAm(a) 



(12) 



where f2,„ is the fractional matter density today. 

The evolution of an additional dark component could 
be described using phenomenological models for macro- 
scopic variables, such as the dark energy equation of state 
and sound speed [tJ: however this leads to associated 
theoretical priors on the growth history [t^. Here we 
make no assumptions, apriori, about the origin of the 
modifications in Models I and II, whether gravitational 
or from some sort of inhomogeneous matter [76j . 



C. Parameterizing the growth history 

The constraints on the phenomenological model pa- 
rameters give a simple indication of whether deviations 
from GR, described by the modified Poisson equation and 
inequality between (j) and -0, are allowed by the data. It's 
instructive, however, to also consider how these models 
affect the growth of structure, using parameters that are 
more model independent. Here we consider four such, 
gauge-invariant parameters: 

• The total matter power spectrum today 



P{k) = 



(13) 



where I^tot = (PcAc -I- Pb^b)/ (Pc + Pb) is the gauge 
invariant, rest frame, over-density for CDM and 
baryons combined. 

• The rest-frame CDM over-density growth factor 
relative to today, go, 



go{k,a) 



A,(fc,a) 



Ac(fc,a = 1) 



(14) 



• The growth rate of rest-frame CDM over-densities, 
parameterized by the effective power law index /, 



(11) 



/(fc,a) 



din Ac(fc, a) 
din a 



(15) 



In GR / = 1 in the CDM era, and / < 1 at late 
times in response to accelerated expansion. 
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• The change in the growth rate, relative to the 
change in the matter density 



7(fc,a) 



ln/(fc,a) 
lnfi,„(a) 



(16) 



7 takes a well-constrained value « 0.55 at late times 
in ACDM [t^]- Because of this it, in particular, 
has been widely discussed as a useful parameter 
with which to search for deviations from a standard 
growth history [7S 



their luminosities [S^, and Ag = 1.7h~^Mpc fits numer- 
ical simulations of the galaxy power spectrum [8l| . 

The transfer functions for ISW, galaxy and weak lens- 
ing sources are 



Ti = e 




(19) 
(20) 
(21) 



III. DESCRIPTION OF DATASETS USED 

To constrain the background expansion history, we use 
the combined analysis of the 'Union' set of supernovae 
from the Supernovae Cosmology Project js^ ]. and ad- 
ditional supernovae sample of [55|, the joint analysis of 
BAO acoustic signatures in the Two-Degree Field (2dF) 
and DR7 Sloan Digital Sky (SDSS) surveys H^, and the 
WMAP 7-year CMB temperature and polarization data 

To constrain the growth of structure, we consider the 
large scale ISW measurements from WMAP7, the mat- 
ter power spectrum from the SDSS DR7 release [111 , the 
2MASS and SDSS Luminous Red Galaxy (LRG) galaxy 
auto-correlations and ISW-galaxy cross-correlations se- 
lected by Ho. et. al. and the publicly available data 
from the COSMOS weak lensing survey [63. We provide 
more detail on how these datasets are employed here. 

Under the Limber approximation, the 2D angular 
power spectrum for the correlation between two fields, 
X and y, is 



XY 



^Wx{x)WY{x)Tx{ki,x)TY{kux)^l(ki), 
X 



(17) 



where X, Y— I,g and k for ISW, galaxy and lensing con- 
vergence fields respectively, Wx is the window function 
associated with the field X, and Tx is the transfer func- 
tion, with ki ^ {1 + l/2)/x- 

For the ISW and galaxy correlations, we use data 
compiled by Ho et al.H for 2MASS and SDSS LRG 
surveys, with each broken into 4 and 2 redshift bins, 
2MASS0-3 and LRGO-1, respectively. These are cross 
correlated with WMAP5 CMB temperature data. Fol- 
lowing [5^, we model the non-linear corrections and 
redshift-dependcnt bias for each galaxy bin (for the cor- 
relations in the Ho et al. samples only) using a Q-model 
for the matter power spectrum [8l| . 



Pobsik) 



.1 



1 



Agk 



Punik), (18) 



brei{z) is a fixed, unit normalized function describing the 
redshift variation of the galaxy bias in each bin based on 



where Treion is the reionization optical depth, X is the 
transfer function of X, normalized so that X'^{k,x) = 
X^(A:, x)A^(A:) where A|j(fc) is the dimensionless pri- 
mordial spectrum of curvature fluctuations. 

The galaxy and lensing window functions are depen- 
dent on the distribution of galaxy number density in each 
redshift bin i for the relevant survey, n^{x)i 



wlix) = K,i{xW{x) 



x) 



(22) 
(23) 



where 71* is normalized such that, dxn^ix) = 1- The 
bias- weighted mean redshifts (z) — J dzbrei{z)n{z) for 
the 2MASS0-3 and LRGO-1 samples are 0.06, 0.07, 0.10, 
0.12, 0.31 and 0.53 respectively. For ISW, Wi{x) = 1. 

Following [5^, as summarized in Appendix |B1 we 
treat bnorm and Qg as free parameters that are ana- 
lytically marginalized over [83| using the galaxy auto- 
correlations. This optimal value of the bias is then used 
to normalize the ISW-galaxy cross-correlation spectrum. 



We use the publicly available analysis of the COSMOS 
weak lensing survey data [60l | , given in terms of the pro- 
jected power spectrum of 2D shear correlations. 



1 

2^ 



Co Ad) 



dl ICr-hA^d)- 



(24) 
(25) 



The data is broken up into 3 redshift bins, z < 1, 
1 < z < 1.4 and 1.4 < z < 3. To obtain the weak 
lensing convergence power spectrum, C;'*", we apply the 
Smith et al. |84l | non-linear correction to the linear power 
spectrum. We use the weak lensing likelihood code writ- 
ten by Julien Lesgourgues [ssl ] and follow their modeling 
of the systematic errors in the weak lensing data. Three 
systematic errors parameters are explicitly marginalized 
over: A accounts for a 6% uncertainty in the overall cal- 
ibration of the shear measurement based on simulated 
HST images, it has a Gaussian prior centered on 1 with 
a A = 0.06; B accounts for a 5% relative calibration un- 
certainty between the shear measured in the high and 
low redshift bins due to underestimation of shear in faint 
or small galaxies, it has a one-sided Gaussian prior with 
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FIG. 1: The efTect of changing Q [left panels] and R [right 
panels] in Model I, defined in ([8]). Shown are the CDM density 
fluctuation, Ac, its growth rate, /, and ISW and weak lensing 
transfer functions, Ti and Tk, today, relative to their values 
for ACDM. Four values of Qq and i?o are shown for a time- 
varying but scale-independent modification with s = 3 and 
fcc = oo, and all other cosmological parameters fixed. On 
subhorizon scales, boosting the gravitational potential (Q > 
1) boosts the growth of dark matter density perturbations 
and the lensing potential. If Q is increased, on subhorizon 
scales the boost in the gravitational potential opposes the 
late-time decay induced by accelerated expansion. The ISW 
signal can be progressively reduced, canceled out (for Q a 
little larger than unity) and turn negative as Q increases. 
The modification has a larger and opposite effect on CDM 
growth on horizon scales, where the time evolution of the 
modification plays a key role. Increasing Q causes growth in 
Ac to slow, or even decay at late times on the largest scales 
(with / = dlnAc/dlna < 0). The effect of increasing R 
is quantitatively degenerate for subhorizon CDM growth and 
qualitatively similar otherwise, but with a slightly reduced 
amplitude. 



FIG. 2: The effect of changing A^h = 0((S'7?),±5 x 10"^ 
in Model II defined in for two different scale depen- 

dences, Tia; = O.StIs [left panels] and n^; = 0.75ns [right pan- 
els]. For this illustration we fix Sa, = 1 and include a cut 
off, — O.OlMpc""^, and all other cosmological parameters 
fixed. Shown are the CDM density fiuctuation, Ac, its growth 
rate, /, the ISW transfer function, Tj, and equivalent value 
of (3, today, relative to their values for ACDM. A positive 
Aj,H is equivalent to (5 > 1, suppressing cold dark matter 
growth, that can give rise to / < on large scales, and boost- 
ing the ISW component. For n^^ = 0.5ns the modification 
has roughly the same spectral evolution as the CDM pertur- 
bations, leading to a scale independent effect on large scales. 
Increasing the spectral tilt, n^, exacerbates the impact of the 
modification for scales just below the horizon. 



IV. OBSERVATIONAL IMPACT OF A 
MODIFIED GROWTH HISTORY 

We can combine the CDM fluid equations, ([2]) and ([3]) 
for w = = 0, to give a second order evolution equation 
for Ac that makes no assumptions about the underlying 
gravitational theory, 

Ac "f H A, + M = 



(Jb = 0.05 and i? > 1; and a 10% intrusion of low-rcdsliift 
galaxies into the high redshift bins due to catastrophic 
photometric redshift errors is modeled by a third parame- 
ter, (7, using a a one-sided Gaussian prior with ac = 0.10 
and C > 1. The correlation functions Ci^2 are multiplied 
by (A/B)'^, and {AB)'^/C respectively in the low, 
medium and high redshift bins. 



with 



-3[^ + nij + nh-:^\{(j) + n^) + n^p\ (26) 



/i = fc2 + 12nGa^J2(P^ + (27) 

i 

On sub-horizon scales (1261) reduces to the commonly 
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FIG. 3: The effect of the scale-independent and time-evolving 
modifications in figure [T] on the CMB [upper panel] , lensing 
correlation [middle panel] and ISW-galaxy cross-correlation 
[lower panel]. For the lensing and ISW-galaxy correlations 
example redshift bins are shown. With increasing Q, the lens- 
ing potential, and galaxy number density are boosted while, 
by contrast, the late-time ISW is suppressed. Both positive 
(Q < 1) and negative (Q > f) ISW amplitudes give rise 
to a boost in the CMB power spectrum at large scales. In- 
creasing Q reduces the ISW-galaxy cross-correlation, and can 
lead to anti-correlation at low redshifts (this is relevant for 
the 2MASS samples, but not the higher redshift LRG sam- 
ple shown here). Increasing R has a qualitatively degenerate 
effect on the spectra. 



considered equation for CDM growth, 

Ac + nAc + k^iJj w 0. (28) 

In the absence of anisotropic shear stresses, k'^ip « 
—AirGa^QRA. Individually increasing Q oi R would 
have degenerate effects on Ac, countering the suppres- 
sion in growth of CDM over-densities brought on by cos- 
mic acceleration, increasing the lensing shear correlations 



FIG. 4: The effect on the CMB temperature power spectrum 
of a time-independent modification to the Poisson equation 
[top panel] and the relationship between the two potentials, 
(fi and V' [lower panel] . The principle effect of a larger Q is to 
suppress the early ISW effect on scales comparable with the 
horizon scale at last scattering. A larger R generates a late- 
time ISW amplitude even in the matter dominated era, by 
introducing a multiplicative difference between ip and tp and 
boosts the temperature dipole anisotropy, raising the acoustic 
peaks. In contrast to the time varying evolution in figure [S] 
constant Q and R have opposing, rather than degenerate, 
effects on the CMB. 



and reducing the ISW component of the CMB tempera- 
ture perturbation. As Q or i? increases, the correlation 
between ISW and the galaxy distribution decreases and 
can give rise to an anti-correlation at late times. 

On super-horizon scales, the CDM evolution is deter- 
mined by the time derivatives of the potentials, 

Ac + 2HAc - 3 + H{2^ + V') + (S'H + ■H^)i?j w 0. (29) 

On these scales, a modification that boosts the potential 
can lead to a decline in the growth rate of Ac and the 
lensing potential, and an associated boosting of the ISW 
transfer function. 

Figure [T] illustrates the effect of the modifications, Q 
and R, on sub- and super-horizon scales for Models I. 
One finds a marked contrast of the responses of Ac and 
the transfer functions on large and small scales consistent 
with the discussion above. On small scales boosting the 
potential promotes CDM growth, boosting the lensing 
transfer function, and suppressing the ISW effect. Large 
scales show the opposite behavior, with CDM growth 
suppressed by a boost in the potential. Note that on 



7 



horizon scales suppression of CDM growth can give rise 
to a negative growth rate, / < 0, for which the 7 pa- 
rameter is undefined. The opposing behaviors on large 
and small scales suggest that measuring the growth his- 
tory on all scales could reveal distinctive signatures of a 
modified growth history. 

Figure [5] illustrates how having the additive modifica- 
tion to the Poisson equation, in Model II, modifies the 
growth history on large scales. Qualitatively the effect 
is similar to that for Model I, dependent on the sign of 
Q. AxH > gives Q > 0, suppressing CDM growth and 
boosting the ISW effect. A key difference between Mod- 
els I and II is the spectral dependence in Model II, pa- 
rameterized by n^;. A spectral index, sa ng/2, matches 
the scale behavior of CDM for a scale invariant primor- 
dial spectrum and hence just acts to give a scale inde- 
pendent boost to growth on large scales. Introducing a 
stronger scale dependence accentuates the impact of the 
modification on sub-horizon scales (for scales above the 
cut-off scale kx). This can lead to over-densities being 
transformed into voids for scales just below the horizon, 
giving rise to a sharp asymptote and change in sign in Q, 
as Ac transitions from positive to negative. This is dis- 
tinctly different to the monotonically varying Q in Model 
I. 

In figure [3] we illustrate the effect of a time-evolving 
modification on the CMB, Icnsing and ISW-galaxy cor- 
relations. Since the CMB power spectrum is effectively 
related to the magnitude of the ISW signal, the spectrum 
docs not have a monotonic response to Q and R] there is 
a minimum ISW contribution for Q or R a little greater 
than unity, when the suppression in gravitational poten- 
tial growth rate due to cosmic acceleration is most closely 
countered by the boosted from the modificatioii, an ef- 
fect seen in previous modified gravity analyses (FTI . [67| . 
For larger values of Q or i?, a negative late-time ISW 
signal gives rise to a boost in the large scale CMB power 
spectrum. The lensing power spectrum is similarly not 
monotonically responsive to Q and R, having a minimum 
for i? = — 1 ; however this degeneracy is removed through 
our theoretical prior R> —1. 

The CMB power spectrum alone cannot distinguish 
between boosting > 1) or suppression {Q,R < 1) 

of the gravitational potential since both boost the large 
scale power. ISW-galaxy correlations vary monotoni- 
cally with Q and i?, however, and break this degener- 
acy. Boosting the gravitational potential promotes the 
formation of large scale structure while suppressing the 
ISW signal, hence as Q or i? increase there is a mono- 
tonic suppression of the cross-correlation that can lead 
to anti-correlation at low redshift. 

The late-time modification to the growth history in 
this analysis leaves a degeneracy between Q and R un- 
broken (the data is principally sensitive to the sum of 
the potentials on smaller scales, or its rate of change on 
larger scales, sensitive toQ{l + i?)/2). An increase in Q 
can be countered by a decrease in R. In order to break 
this degeneracy additional correlations would have to be 



TABLE I: ID marginalized 95% confidence limits for modi- 
fied gravity parameters Q and R, and principle degeneracy 
direction Q — R, for time- and scale-independent deviations 

from GR in Model I. 

Data analyzed Q R Q — R 

All [0.97, 1.01] [ 0.99, 1.02] [-0.032, 0.009] 

No WL [0.96, 1.01] [ 0.98, 1.02] [-0.034, 0.008] 

No WL, no ISW-galaxy [0.98, 1.04] [ 0.99, 1.03] [-0.025, 0.014] 



employed, such as using galaxy-peculiar velocity correla- 
tions or galaxy-lensing correlations, and considering ra- 
tios of observables [1, • 

In figure |4] we show how time-independent Q or R have 
opposing, rather than complementary, effects on the first 
CMB acoustic peak. In the standard GR scenario during 
the matter dominated era, the gravitational potential is 
constant and no ISW signal is generated. An early ISW 
effect is present around recombination however, because 
the universe is not purely matter dominated causing the 
gravitational potential to decay. This decay boosts the 
CMB power spectrum on scales of order the horizon size 
at last scattering. A constant Q > 1 counters the early 
ISW and suppresses the CMB temperature correlations 
below a degree scale. Because it is constant, and parity 
between and ip is maintained at late times, it has an 
insignificant effect on the late time ISW. Increasing R 
has two effects: firstly it introduces a constant disparity 
between the growth rate of the two potentials which leads 
to a non-zero ISW signal even during the matter domi- 
nated era, and secondly it increases the dipole anisotropy 
in the CMB, boosting the acoustic peaks. In contrast to 
the late-time modifications to the growth history, any in- 
crease in a time- independent Q could be countered by 
an increase in R, to minimize the impact on the well- 
measured height of the first peak. We willfind that, con- 
sistent with this, the data effectively measures the com- 
bination Q — R. 



V. MCMC ANALYSIS 

We perform Monte Carlo Markov Chain (MCMC) 
analyses using a version of the CosmoMC code [8j| modi- 
fied to include general growth evolution models described 
in section |lll In scalar tensor theories of gravity, such 
as f{R), one can construct the modification to GR that 
exactly matches a chosen background expansion history 
[SmSgf . but gives rise to modifications of the growth 
of perturbations (Tlj . With this in mind we assume a 
standard ACDM expansion history and focus on the con- 
straints on the modified growth history. 

We impose flat priors on six core cosmological pa- 
rameters: the baryon matter density, ^Ibh^, the CDM 
matter density, ric^i^, the angular size of the last scat- 
tering horizon, 0, the reionization optical depth, r^ei, 
the primordial scalar power spectrum index and ampli- 
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TABLE II: ID marginalized 95% confidence limits on the modified gravity parameters Qo, i?o and Qo 
principle degeneracy direction, well described hy Q{1 + R) /2. 
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FIG. 5: 68% and 95% confidence limits for the time- and 
scale-independent modifications in Model \ (Q = Qo and 
i? = i?o, s = and kc — oo). Constraints are shown for all 
data combined (forefront, yellow/orange contours), excluding 
the lensing correlation (red, dot-dashed) and excluding both 
lensing and galaxy-ISW correlations (blue, dashed). The ef- 
fect of a constant modification on correlations around the first 
acoustic CMB peak leads to tight constraints, limited to de- 
viations < 3% from GR at the 95% confirndence level for 
all data combined. Driven by their opposing effects on the 
well-measured first peak, the principle degeneracy direction 
is well-described by constant {Q — R). 



FIG. 6: As in figure[5]but giving the 68% and 95% confidence 
limits for allowing scale-independent, but time-evolving mod- 
ifications in Model I, that become important at late-times 
(0 < s < 3, fcc = oo). The principle degeneracy direction 
is well-described by constant Qo{l + Ro)/2. Constraints on 
the time-evoving modification are primarily driven by their 
effect on the cosmic-variance limited, late-time ISW effect in 
the CMB. As such, they are improved by the inclusion of 
the galaxy-ISW cross-correlation data, but remain far weaker 
than for the time- independent modification shown in figure [5] 



tude, Us and InlO^^As respectively. Here h = _ffo/100 
kms~^Mpc~^ where Hq is Hubble's constant. Following 
[ssj , we marginalize over absolute calibration uncertain- 
ties in the 3 weak lensing bins using 3 nuisance parame- 
ters A, B and C as discussed in lllll In each scenario, we 
run 8 independent chains until they satisfy the Gelman- 
Rubin convergence criteria [ooj . 

For Model I, we consider both a timc-dcpcndent and 



time-independent scenario. For the time-independent 
scenario (s = 0) we assume scale invariance (fcc oo) and 
impose flat priors on < Qo.oo < 10 and —1 < i?o,oo < 
10. Table U and figure [5] summarize the marginalized ID 
and 2D constraints on Q and R and the principle de- 
generacy direction. The well measured CMB power spec- 
trum on sub-degree scales provides the dominant con- 
straint, limiting the effect of time-independent modifica- 
tions to Newton's constant , Q-l< 3% those from stan- 
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dard gravity at the 95% confidence level. This is compa- 
rable with constraints on time- independent modifications 
to the Newton's constant from bounds on the expansion 
history with the CMB [U . The galaxy-ISW and lensing 
correlations only provide limited improvements. As dis- 
cussed in section HVl R and Q have opposing effects on 
the first peak and the degeneracy direction of the con- 
straints is consistent with this, being well described by 
Qo — Rq consistent with zero, and both Qo and Rq con- 
sistent with unity at the 95% confidence level. 

We also investigate time-dependent scenarios for 
Model I, in which the modification characteristically 
arises at late times. Here we utilize the degeneracy di- 
rection Q{1 + R) /2 to more efficiently search the param- 
eter space, imposing flat priors on < Qo.oo < 10 and 
< Qo,oo(l + ^o,oo)/2 < 10, and rejecting combinations 
with R < —1. To ensure the choice of prior did not 
unduly affect the constraints we also ran MCMC chains 
using flat priors on Qo.oo ^ [0, 10] and i?o,oo ^ [—1, 10] 
and found no significant change in the results. We con- 
sider both scale independent modifications, with fcc = oo 
(for which Qoo and i?oo arc redundant), and scenarios in 
which large and small scales are able to vary indepen- 
dently. 

Table summarizes the ID marginalized constraints 
on late-time evolving modification for which we marginal- 
ize over < s < 3. For scale- independent modifications 
(fcc = oo), one finds that the constraints on the principle 
measured combination, (3o(l+-Ro)/2, arc only marginally 
affected by omitting the weak lensing data but are more 
significantly relaxed if the galaxy-ISW correlations are 
not included in the analysis. This presence of this prin- 
ciple mode is demonstrated by the insensitivity of the 
constraint to whether Qo or Rq alone, or both, are al- 
lowed to vary, and is clearly shown in Figure |6l 

For scale-dependent modifications we consider a tran- 
sition scale of kc = lO^^Mpc^^. The constraints on 
the large scale modifications arc significantly more re- 
laxed than on the smaller scales, with the dominant con- 
straints coming from the CMB ISW. As shown in Fig- 
ure [3 the data effectively measures Qo(l + ^o)/2 and 
Qoo(l + i?oo)/2 at large and small scales, respectively. 

The change in the maximum likelihoods of the sce- 
narios considered in Model I reflect the degeneracies in- 
herent in the constraints. For all data combined, the 
flducial ACDM scenario has bcst-flt likelihood, C, given 
by — 21n£ = 4070.3. The scale independent scenario in 
Model I has a best-flt —2 ln£ = 4067.7; though 3 param- 
eters are varied, Qq, Rq and s, only one measured degree 
of freedom is really introduced, effectively Qo{l + Ro)/2, 
while s is unconstrained. The current data does not in- 
clude sufficient information to constrain any scale depen- 
dence in the modification, reflected in no improvement in 
the best fit in comparison to the scale independent model. 

For the scale-independent and dependent models in 
figure [71 we calculate the matter power spectrum to- 
day in 40 logarithmically spaced bins in k from 10~^ — 
10~'^h/Mpc, and the growth parameters /, go and 7 in 30 



bins of width Az = 0.1 between z = and 3. These allow 
us to find the ID marginalized distributions for the power 
spectrum and growth parameter values in each bin. 

In figure [5] we show the ID marginalized constraints 
on the matter power spectrum for scale independent and 
dependent scenarios. Even in the scenario in which the 
small scales are allowed to vary separately from the large 
scales, there is only a small broadening in the amplitude 
of the spectrum in comparison to that from ACDM. By 
contrast, the matter power spectrum on large scales can 
boosted or suppressed relative to ACDM at the 95% con- 
fidence level, with the 68% confidence levels driven by the 
low, large scale CMB power. 

In figure O we show the constraints on the general 
growth parameters /, go and 7 for two comoving scales, 
k — 0.005/Afpc and 0.05/Mpc. The suppression in large 
scale growth seen in the power spectrum today is also 
seen in a lower growth rate (and higher 7 and go) at 
low redshifts. 7 on large scales can take a dramati- 
cally broader range of values than those predicted by 
ACDM. Interestingly, on the largest scales the growth 
rate, / = dlnAc/dlna, can become negative, leaving 
7 = ln//lnf2m undefined. Though 7 is a good parame- 
terization for indicating a modification to gravity on sub- 
horizon scales, many growth histories consistent with the 
current data cannot be described by 7 for scales just be- 
low, or at, the horizon scale. 

For Model II, we focus on how modifications ef- 
fect large scales, considering a cut-off scale of k^ — 
10~^Mpc~^. The modification principally effects modes 
that have entered the horizon after radiation-matter 
equality. We impose flat priors on Axh ^ [~10~"^, 10"'^], 
and marginalize over < s^, < 2 and < < 3. With 
all data combined, we find ID marginalized constraints 
on A^H = [-1.5 X 10-"^, 3.7 x 10"^] at the 95% c.l. There 
is not a significant improvement in the maximum likeli- 
hood in comparison to ACDM, with — 21n£(Model II) 
=4069.8, reflecting the fact that the constraints are en- 
tirely consistent with GR. 

In figure [TOl we show the effect modifications in Model 
II have on the matter power spectrum and the growth 
parameters /, go a-nd 7. The data allows both suppres- 
sion and growth of large scale power in comparison to 
ACDM and with a greater amplitude than the parame- 
terization in Model I. Scenarios with a late-time decay 
of inhomogeneities, with / < (and undefined 7) at 
k = 0.005/il/pc are consistent with the data at a 95% 
confidence level at low redshifts. 



VI. DISCUSSION 

In this paper, we have presented constraints on the 
growth history of cosmic structure with current cosmo- 
logical data in light of phenomenological modifications 
to the Poisson equation and the relationship between 
the two gravitational potentials, respectively described 
by two scale- and time-dependent functions. Q and R. 
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FIG. 7; Impact of allowing different modifications to the growth history on large and small scales for Model I. The figure 
shows 68% and 95% confidence limits on large [left panel] and small scales [right panel] using all data combined. In each figure, 
constraints are shown for the scale-independent scenario (forefront, red contours), and the scale-dependent scenario in which 
distinct deviations from GR are allowed on large and small scales {Qo, Ro,Qoo and Roa are all varied with cross-over scale 
kc = W~^Mpc~^ and < s < 3) (blue, dashed). 



We consider two parameterizations for these modifica- 
tions. The first, Model I, considers monotonically vary- 
ing Q and R in both time and space. The second. Model 
II, leaves the relationship between the two potentials un- 
changed (i? = 1) but allows an additive adjustment to 
the Poisson equation, that could Model either a modifi- 
cation to how gravity and matter are related, or equally 
the presence of an additional inhomogeneous, dark com- 
ponent; from a phenomenological standpoint each inter- 
pretation is equivalent, since the effect on structure for- 
mation is only detectable indirectly through its impact 
on the gravitational potential. 

We find that time-independent modifications in Model 
I, through their impact on the well-measured CMB acous- 
tic peaks, are tightly constrained by the data. In terms of 
the effect on the Poisson equation the constraints trans- 
late into the effective Newton's constant being within 3% 
of its standard value in GR. 

An evolving modification that becomes important at 
late times, as one might expect if it is associated with 
the onset of cosmic acceleration, is, by contrast, far more 
weakly constrained; an effective Newton's constant in 
Poisson's equation today can multiple times its value 
in GR, rather than a few percent variation, and remain 
consistent with the data. The reason for this is that 
the modifications have their biggest effect on the large 
scales (fc > kf,q) that, at present, principally constrained 
through the cosmic-variance limited ISW and a small 
number of data from ISW-galaxy cross-correlations. The 



inclusion of a wider variety of improved precision, cross- 
correlations on these larger scales, as will come from up- 
coming wide-scale surveys, should greatly improve the 
constraints on these models. To precisely compare ob- 
servations to predictions on these scales will involve con- 
sidering the effect of gauge issues in relating theory to 
observables [o^] and the impact of scale-dependence and 
nonlinearity of the bias on large scales [9J, , both of 
which could obscure the effect of an underlying modifi- 
cation to gravity or clustering dark energy. 

The modified growth histories considered here can have 
markedly different behaviors on large and small scales; a 
boosting of CDM growth on small scales can be tied to 
a suppression on larger scales. This opens up the ques- 
tion of using sub-horizon and horizon scale observations 
separately to directly probe for this distinctive signature. 

There is a degeneracy present in constraints from the 
current data, well described by Q{1 + i?)/2, that pre- 
vents one from establishing whether modifications to the 
growth history arise from Q ^ 1 or R 1. If we hope to 
be able to use observational constraints to phcnomeno- 
logically ascertain what is driving any modification to the 
growth history then breaking this degeneracy is critical. 

The application of tailored estimators built from cross- 
correlation of galaxy, weak lensing and velocity fluctua- 
tions, could directly confront this issue. As such, the 
combination of deep and wide, imaging and spectroscopic 
large scale structure surveys, coming online in the next 
few years, provides the exciting prospect of being able to 
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FIG. 8: ID marginalized 68% and 95% confidence limits on 
the total matter power spectrum today for ACDM (black con- 
tour lines) and the two scenarios using Model I from figure [7| 
(filled red contours). The upper panel shows the constraints 
for the scale-independent modification while the lower panel 
shows constraints if both large and small scale growth histo- 
ries are modified independently. The data allows a broader 
range of growth histories, tending to suppress power on large 
scales relative to the ACDM scenario, constrained by the ISW. 
By comparison, the matter power spectrum on smaller scales 
is only marginally effected by the modified growth history. 
This is even the case when the subhorizon scales are modified 
independently of the superhorizon scales. 



direct test the cosmic growth history, and the relation- 
ship between matter and gravity on cosmic scales. 
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perturbations h and rj, rather than the Newtonian gauge. 

Using the notation of Ma and Bertschinger [6l| the 
synchronous and Newtonian gauge variables are related 
by 



1] — Ha, ^ a + Ha 



where 



k'a 



h 
2 



37). 



(Al) 



(A2) 



The synchronous (s) and conformal Newtonian (c) 
gauge matter perturbations are related by 



ak^ 



(A3) 
(A4) 
(A5) 



GAMB evolves the metric perturbation ry and the mat- 
ter perturbations for each species according to the syn- 
chronous gauge fluid equations. To this end, it evaluates 
the following variables at each time step 



crcAMB — ka = 



ZCAMB 



h 
2k 



?7 - ^ 
H 



C^CAAIB 



(A6) 
(A7) 



where, for the modified gravity model, (j) now comes from 
the new Poisson equation (j6|), and the evolution equation 
for 77 is given by the time derivative of ([5]) coupled with 
dZl), Q and dSl), 



JQ , 



(1 + w 
(i? - \)QH\ A, 



where 
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(AS) 
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Appendix B: Galaxy auto-correlation analytical 
marginalization 



Appendix A: Incorporating Modifications to Gravity 
in the CAMB code 

We summarize here how the modified gravity equations 
discussed in section III Bl are incorporated into GAMB 
[63| . GAMB uses the synchronous gauge, with metric 



This section summarizes the analytical marginalization 
approach for 'nuisance' bias and non-linear fitting pa- 
rameters, h and Qg, used for analyzing the ISW-galaxy 
cross-correlation data, following the approach of [53 • 

Consider d = {d* = Cgg(Z)}, the vector containing the 
galaxy auto-correlation data, with each redshift bin, i. 
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FIG. 9: ID marginalized 68% and 95% confidence limits for Model I for the growth parameters /, go and 7 at two scales 
fc = 5 X 10~^i\fpc~^ and k — 0.05Mpc~^ (left and right panels in each figure respectively). ACDM (red contours) is compared 
to Model I with < s < 3 (blue contours). [Left figure] Constraints for the scale- independent modification, with fcc = 00. 
[Right] Constraints for the scale-dependent modification varying {Qo, Ro,Qoo, Roo} with a transition scale kc = Mpc~^ . 
The modified gravity models can have a markedly broader range of growth rates than ACDM, especially on large scales. 



and C is the associated covariance matrix. To analyti- 
cally marginalize to find bnorm and Q in (|18p we construct 
two theoretical auto-correlation sets 



^w,{xrT,{k,xrAUk){Bi) 
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The values of bnorm and Q that minimize the x^ fit to 
the data are evaluated using the matrix M and vector v, 



with xlq = f^C^iq etc. Then 
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with an overall fit of 



= Xld - v^M-V. 
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